Fluctuation-dissipation theorem and flux noise 
in overdamped Josephson junction arrays 
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The form of the fluctuation-dissipation theorem for a resistively shunted Josephson junction 
array is derived with the help of the method which explicitely takes into account screening effects. 
This result is used to express the flux noise power spectrum in terms of frequency dependent sheet 
impedance of the array. The relation between noise amplitude and parameters of the detection coil 
is analysed for the simplest case of a single-loop coil. 
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I. INTRODUCTION 

The two basic experimental methods used for con- 
tactless investigation of finite frequency properties of 
two-dimensional superconducting systems (such as thin 



Josephson junction arrays |6|-|9|] and wire net- 
10|Jllf|) are two-coil mutual-inductance technique 



1 1 1 and flux noise power spectrum analysis 



films 
works 

The first of them is based on measurement of the voltage 
induced in the detection coil by the currents flowing in 
the sample under the action of ac electric field produced 
by the current in the other (driving) coil. For the given 
geometry of the coils the measured signal can be used 
to extract JIIJ the complex frequency dependent sheet 
impedance Z a (w) of the sample on the assumption that 
for wave-lengths larger than the characteristic dimensions 
of the detection coil Z n (cj) is not wave-length dependent. 

In the case of the flux noise spectrum analysis the ap- 
proaches to interpretation of the experimental data are 
much more varied. The theoretical predictions of the flux 
noise spectrum used for comparison with experimental 
data are found by relating it with Jz|JT^] or (no less of- 
ten) by replacing it by p 13,1^] a correlation function 
describing the vortex distribution and, naturally, turn 
out to be dependent on the particular choice of assump- 
tions concerning the form of this distribution. Numerical 
simulations also demonstrate a clear tendency towards 
studying the vortex number noise |jn|[lj|,[l6) rather then 
the flux noise. The only attempt to achieve a description 
of the flux noise power spectrum in terms of the sheet 
impedance of the sample taking into account the actual 
geometry of the detection coil has been undertaken by 
Kim and Minnhagen jl7|. However, this calculation is 
also based on expressing all quantities in terms of the 
vortex gas correlation functions and, therefore, can be 
trusted only in a limited range of parameters. 

In the present work we argue that in the case of a 
resistively shunted Josephson junction array the general 
expression for the flux noise spectrum can be found with- 
out artificial decomposition of all fluctuations into the 
vortex part (which is usually assumed to be responsi- 
ble for the flux noise) and the remaining so-called " spin- 



wave" part (which is traditionally neglected). Although 
in semi-phenomenological treatment p2]- [l4| , |T8| ] of two- 
dimensional superconductors such decomposition seems 
to be inevitable, the case of an overdamped Josephson 
junction array allows for application of a more universal 
approach to calculation of the flux noise power spectrum. 
It is based on the direct relation (discussed in Sec. 2) of 
the flux noise with current correlations, which, on the 
other hand, can be expressed in terms of the complex 
frequency dependent sheet impedance with the help of 
the fluctuation-dissipation theorem. 

The additional advantage of such approach is that it 
allows to include into consideration in systematic way the 
mutual influence between magnetic field fluctuations and 
current fluctuations (the screening effects), which insofar 
has been neglected in theoretical works |l2]-|l7| devoted 
to flux noise spectrum analysis. The form of the Hamilto- 
nian, which should be used for the description of a resis- 
tively shunted array in presence of self-induced magnetic 
fields, is discussed in Sec. 3, and the corresponding dy- 
namic equations in Sec. 4. 

The explicit form of the fluctuation-dissipation theo- 
rem for resistively shunted Josephson junction array is 
derived in Sec. 5. It shows that the current correlations 
in the array are determined by the response of the cur- 
rent to the external electric field and not directly by 
the sheet impedance of the array (which is defined as 
a response to the total electric field) The nature of the 
expression for currents correlation function, related with 
the peculiarities of the two-dimensional geometry, allows 
to expect the same expression to be applicable for arbi- 
trary two-dimensional systems in which capacitive effects 
can be neglected. 

Our main result, the relation between the flux noise 
power spectrum and the frequency dependent sheet 
impedance of a two-dimensional superconductor, is pre- 
sented in Sec. 6, which includes also the discusssion of 
the noise spectrum dependence on the parameters of the 
detection coil and comparison of the results with those 
of other authors. 
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II. FLUX NOISE AND CURRENT 
CORRELATIONS 

In a flux noise experiment one measures and analyses 
the time dependence of a voltage created in a detection 
coil by fluctuations of currents in some conducting (or 
superconducting) object. This voltage is determined by 
the time derivative of the magnetic flux penetrating the 
coil, and the value of the flux can be expressed in terms 
of the current density distribution j(r) [r = (x%, X2, 2:3)] 
inside the object with the help of the Biot-Savart's law, 
which in the Coulomb gauge (divA = 0) can be written 
as 



AA(r) = - M oj*(r) 



(1) 



where A(r) is the vector potential defining the dis- 
tribution of the magnetic field (magnetic induction) 
B(r) = rot A(r) created by j(r), 



A = 



d 2 d 2 d 2 



(2) 



dx\ dx\ dx\ 

is the three-dimensional Laplacian and 

j t (r)=j(r)-A- 1 g raddivj(r) (3) 

is the transverse part of j(r). Magnetic fields produced 
by the longitudinal part of j(r) cancel each other. 

In the case of a system which can be considered as ef- 
fectively two-dimensional and situated (for simplicity) in 
the plane x 3 = the three-dimensional current density 
j(r) is reduced to 



j(r) = i(x)J(x 3 ) 



(4) 



where x = x a (a = 1, 2) is the two-dimensional vector 
defining the position of a point in the plane x 3 = and 
i = i a is the two-dimensional vector describing the two- 
dimensional current density. 

Substitution of Eq. ^ into Eq. (|l|) allows then to find 
that 



A(r) = ^0 



d 2 q 



dq 3 expi(qx + q 3 x 3 ). t 



Mo 
2 



(2tt) 2 J 2tt q 
d 2 q_ exp(iqx-g|x 3 |). 
(2n) 2 



where 



i(q) = J d 2 x exp(— iqx)i(x) 



i £ (q) 
i'(q) > (5) 

(G) 



is the (two-dimensional) Fourier transform of i(x), 
q = |q| and 



i*(q) = i(q)-q(qi(q)) 



(7) 



is the transverse part of i(q), q = q/q being the unit 
vector parallel to q. 



In the simplest case a coil can be approximated by a 
closed circular ring. Integration of A(r) over the perime- 
ter of the ring x\ + x\ = r 2 situated at the distance h 
from the plane x 3 — gives 



$ = /<irA(r) = ^ J ^_F(q)z'(q) 



where 



a, 13 



(8) 



(9) 



is the amplitude of i*(q), £ a p is the unit antisymmetric 
tensor, 



2TrrJ 1 (qr) 
F(q) = exp(-qh) 



(10) 



is the geometrical factor depending on the parameters of 
the coil and Ji(z) is the first order Bessel function. In 
the case when the coil can be considered as consisting 
of N turns separated by the distance b from each other, 
the expression for F(q) should also include an additional 
factor obtained by summation of contributions from dif- 
ferent turns JToj ] : 

2itrJ 1 (qr) 1 - exp(-Nqb) 

^(q) = exp(-qh)- — — . (11) 

q 1 — cxp(— qb) 

The power spectrum of the flux noise is given by the 
flux-flux correlation function 



S(uj)= J dt($(t + t)$(t )) exp(iut) (12) 

and with the help of Eq. (||) can be expressed in terms 
of the current density correlation function 

(i*i*)qa, = J d 2 x J dt [exp(-zqx + iujt) (13) 
x (i*(x +x,£ + t) i*(x ,to))] 



as 



(14) 



III. THE HAMILTONIAN OF A JOSEPHSON 
JUNCTION ARRAY 

When self-induced magnetic field is taken into account 
a square Josephson junction array can be described by 
the Hamiltonian fl9| , po| 



H = -J22 cos (Va^n - A aa ) 

n, a 

+^(VxA)„M nk 1 (Vx A) k 



(15) 
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where <p a is the phase of the order parameter on the n- 
th superconducting island, the variables A na (defined on 
the bonds of the lattice) are determined by the integral 
of the vector potential A(r) over the line connecting the 
geometrical centers of the neighboring superconducting 
islands: 



2e / ,a ( n + e °) 
A na = -r / dr A(r) 

n .)„„ 



(16) 



and a is period of the lattice. 

The first term in Eq. ( p"5| ) describes the Josephson en- 
ergy of the junctions in the array. The coupling constant 
J entering this term is determined by the critical current 
I c of a single junction: 



Eq. (|H)] Vx stands for epa^p, whereas V/3 desig- 
nates the lattice difference, analogous to the one defined 
by Eq. ([La), but shifted in the negative direction: 



V/3X n = X n — X n - 



(24) 



On the other hand variation of Eq. (15) with respect 
to ip a gives the current conservation equation 



(V/) r 







where 



(V-T)n = y^[I n q - I(n-e a )a] 



(25) 



(26) 



J = Trie , 
2e 



(17) 



which is assumed to be the same for all junctions, whereas 
V a y n denotes the difference of y>„ between the neighbor- 
ing sites of the lattice: 



<Pt. 



(18) 



Here ei = (1,0) and e2 = (0, 1). Notice that the combi- 
nation 



(19) 



which enters as the argument of the Josephson energy 
Ej{6) = — J cos 9 is a gauge- invariant quantity. 

The second term in Eq. (|l^) is the energy of the mag- 
netic field 



E mi ^^- f d 3 rB 2 (r) 
2^0 J 



(20) 



expressed in terms of the variables A na . The matrix 
M n t = M(n — k) is usually called the mutual inductance 
matrix li9H2l and 



(V X A) n = ^ tpa^pAna 

a,/3 



(21) 



is the directed sum of the variables A- 3a along the perime- 
ter of a lattice plaquette (the lattice equivalent of rot A) 
and is proportional to the magnetic flux penetrating the 
plaquette. 

Variation of Eq. ( |l5| ) with respect to A na gives the 
equation 



2e 



\ " 



h ^ 

k 



VxM.'fVx A) k 



(22) 



which relates the value of the superconducting current 
through a junction 



Ina = Ic sin(V Q </7„ - A na ) 



(23) 



with the vector potential of the magnetic field induced 
by the presence of the currents in the array. Here [like in 



is the lattice equivalent of divergence. Eq. (|25|) can be 
alternatively obtained by application of the operator V Q 
to Eq. @. Therefore Eq. (g|) and Eq. (||) [both ob- 
tained by variation of Eq. Jig)] are not independent of 
each other. 

The vector potential of the magnetic field created by 
the currents flowing in the array can be chosen purely 
transverse (divA = 0), which in terms of A na corre- 
sponds to 



V a j4 nQ — . 
In that case Eq. ( p2] ) is reduced to 
2c 



(27) 



(28) 



where = ^ P is the two-dimensional lattice 
analog of the Laplacian: 



13 



(29) 



Comparison of Eq. (|28|) with Eq. 
for In - k| > 1 



nk 



2e J 7r^oa| n — k| 



allows to find that 



(30) 



whereas for |n — k| ~ 1 the form of depends on the 
particular shape of superconducting islands |2~jJ] . 

Linearization of Eqs. (^) and their solution allows to 
show that when the magnetic fields of the currents in the 
array are taken into account, the logarithmic interaction 
of vortices becomes screened at so-called magnetic 
field penetration length A, exactly as it happens in super- 
conducting films j22j . When screening is relatively weak 
(that is when A ^> a) , the value of A is given by 



A 



/ioJ \2e 



(31) 
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and does not depend on the shape of superconducting 
islands forming the array [fl9| . 

Instead of considering Hamiltonian (|l5|) as dependent 
on two different types of variables defined on the sites 
(tpn) and on the bonds (A aa ) of the lattice, it is conve- 
nient to use a single variable, namely the gauge invariant 
phase difference 9 na defined by Eq. . In terms of 9 na 
the Hamiltonian (|l5| ) can be rewritten as 



H 



2^ 

n,k 



(Vx«) tt J^(Vxfi) kl 



(32) 



variation of which with respect to 9 na reproduces Eq. 
(E2T) in the form 



J - = ^lE VX M nk(V X 0) k 



(33) 



where the expression for the superconducting current 

Ina = I C Sin Ona (34) 

is naturally consistent with Eq. (P3| ). As previously, the 
current conservation equation (|25j) can be obtained by 
application of the operator V Q to Eq. (|33|). 



IV. DYNAMIC FLUCTUATIONS IN ARRAY OF 
RESISTIVELY SHUNTED JUNCTIONS 

The dynamic description of the same system requires 
to complement the Hamiltonian H by the dissipative 
function W (we assume that the array is overdamped 
and therefore its dynamics is purely relaxational). In the 
case of the array formed by SNS (superconductor - nor- 
mal metal - superconductor) junctions one can describe 
dissipation in terms of the effective resistance shunting 
each junction (so-called RSJ-modcl). This corresponds 
to W{9} of the form 



W 



d 

Ft' 



(35) 



where the effective viscosity 



V 



1 

R 



(36) 



is determined by the value of the shunting resistance R, 
which is assumed to be the same for all junctions. For W 
of the form (|3^) the conservation of energy is achieved 
when the time evolution of the variables 9 na is governed 
by the standard equations of relaxational dynamics: 



dH 
d9 na 



(37) 



On the other hand, Eq. (|37]) can be rewritten in the 
form (|33|), where the expression for the current should be 
replaced by 



I na I c sin 9 na 



n a 

2^Rdt 



(38) 



The time derivative of 9 na being proportional to the volt- 
age, the second term in Eq. ( |38| ) can be easily identified 
as the normal current flowing through the junction. Con- 
sideration of purely relaxational dynamics means that 
we are neglecting capacitive effects and currents have to 
be conserved on each site of the lattice (in other words, 
only transverse current are allowed). This is ensured by 
the form of Eq. (|33| ) , substitution of which into Eq. j2q ) 
automatically leads to its fulfillment for any form of I na . 

In presence of thermal fluctuations the right-hand side 
of Eq. ( |37| ) should be complemented with the random 
force term £ nQ (£): 



d 



dH 



at o9 na 



(no "i" /no 3 



(39) 



the correlations of which are Gaussian and satisfy 

(fna(*)&/j(t')> = 2r)T6 nk 5 af3 6(t - t') , (40) 

where T is the temperature expressed in energy units 
(that is multiplied by the Boltzmann constant fcg). We 
also have included in the right-hand side of Eq. ( |39| ) the 
non-random external force / nQ (to be discussed later). 

In terms of the expression for the current the inclu- 
sion into Eq. ( |39| ) of the random term £ nQ corresponds 
to appearance in Eq. 
contribution to norma 



of the additional (fluctuating) 
current SIr, a : 



I r sin ( 



h d 
2eR dt 



SIna 



(41) 



where SI na = —(2e/h)£ yIla . Note that since Eq. ( |33| ) de- 
scribes the relation between the currents and the mag- 
netic field induced by them, it has to remain fulfilled also 
when fluctuations of currents are taken into account. The 
validity of the current conservation equations (^5|) re- 
mains ensured by the form of the right-hand side of Eq. 

The suggestion to describe the dynamics of a resis- 
tively shunted Josephson junction array by Eqs. ( |4l| ) 
has been put forward by Shenoy p3| , who did not con- 
sider fluctuations of the magnetic field, that is assumed 
9 na = (p n + ea — (fin- In that case substitution of Eqs. ( |4l| ) 
into current conservation equations ( p5| ) leads to dynamic 
equations for ip n with non-local effective viscosity [ p3| . 
Quite remarkably the inclusion into consideration of the 
magnetic field fluctuations leads to simplification of the 
dynamic equations which (in terms of 9 na ) become local. 
The idea that in presence of magnetic field fluctuations a 
resistively shunted Josephson junction array can be de- 
scribed by Eqs. (^), where 9 na is the gauge invariant 
phase difference, has been introduced by Dominguez and 
Jose M. 
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V. THE FLUCTUATION-DISSIPATION 
THEOREM 



L-^VxM^Vx A°) k . 



(49) 



It is well known that when the time evolution of some 
variables {9} is governed by the standard Langevin equa- 
tions of the form (|3^) , the equilibrium (that is calculated 
for /no = 0) correlation function 



Cna,k/3(£ — t ) = (8 na (f)#k/3 (t )) f=Q 

is related with the response function 



<5 7kj8 



(42) 



(43) 



by the fluctuation-dissipation theorem: 

1 9 

G nQ: k^(0 - Gna,k0[—t) = ~ ^ ^C'na.k/SW • (44) 

However, in practical situation one is interested not in 
the response of the gauge invariant phase difference 9 na 
to the (unspecified) conjugate force / na , but rather in 
more readily observable quantities such as conductivity, 
which is the response of a current to application of elec- 
tric field. In situation when electric field E(r) is created 
due to presence of ac magnetic field it is given by 



E(r) = -| A (r) 



(45) 



where A(r) is the vector potential defining the (total) 
magnetic field B(r) = rot A(r). 

In presence of the external magnetic field 
B c (r) = rot A e (r) the expression ( |20| ) describing the 
magnetic field energy should be replaced by the ex- 
pression for the Gibbs free energy 



F mi = ^-J rf 3 r[B 2 (r) - 2B(r)B c (r)] , 



(46) 



variation of which in absence of any other terms gives 
B(r) — B°(r). This leads to appearance in the Hamilto- 
nian (f32|) , describing the array, of the additional term 



AH = ^(V x 8) n ^ l (V x A c ) k 



Here the variables 



na h 



a(n+e a 



drA e (r) , 



(47) 



(48) 



defined on lattice bonds, are related to the vector poten- 
tial A c (r) of the external magnetic field exactly in the 
same way as earlier introduced variables A na are related 
to the total vector potential A(r). 

The form of the correction to the Hamiltonian given 
by Eq. (^) corresponds to the presence in Eq. ( ^9| ) of 
the external force term 



Comparison of Eq. ( [49|) with Eq. ( |33|) shows that (up to 
the factor of 2e/%) f na is related to A^ a exactly in the 
same way as I na is related to 9 na . This allows to conclude 
that the correlation functions of the currents in the array 



cL MP (t-t')^(i na (t)i^(t')) 

and the response function 



r i ( . ,r\ _ S{I aa ) 



(50) 



(51) 



have to satisfy the relation 



h 1 d 



nQ,k/3(*) — ^na,kff( *) ~ ^ j! ~QjpliaMP (t) (52) 



completely analogous to Eq. (Q). Here {A^gj 1 is the 
transverse part of A^a. As can be seen from the right- 
hand side of Eq. ([49|), the longitudinal part of A^ g is 
completely decoupled from fluctuations of 9 na . 

In terms of effective conductivity g^ a k o, defined as the 
response function 



sv. 



(In«}\ 



v°=o 



of the current with respect to external voltage 
h d 



|) can be rewritten as 



Eq. 



Cna,k^(0 — ^[3na,k^W + 9na,k/3 ( — *)] 



(53) 



(54) 



(55) 



For wavevectors small in comparison with a -1 , the 
variables I na can be identified with ai(x), whereas V£ a 
with <zEjj(x), where Ejj(x) is the projection of the ex- 
ternal electric field E c (x) on the plane X3 = (here and 
below we assume that Ejj (q, ui) is transverse) . This allows 
to rewrite Eq. (Kq) as 



(i t i t ) qw = 2TRe[ S eff (q,c)] , 



(56) 



where the Fourier transform g (q,w) of the effective 
conductivity is the coefficient of proportionality in the 
relation 



i(q,w)= ff cff (q, W )EUq, W 



(57) 



One should distinguish between g cff (q, ui) and (also 
momentum and frequency dependent) sheet conductance 
gu (q, ui), which is defined as the coefficient of proportion- 
ality between i(q, ui) and total electric field E||(q, ui): 



i(q,w) = 5 D q,w)E||(q,w) 



(58) 
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The form of the current induced correction to Ejj (q, to) 
can be easily found with the help of Eq. (||) and Eq. (45), 
which lead to 



E||(q,w)=E|(q,w) 



• Mo . / \ 
i(q,o;) 

2q 



(59) 



Substitution of Eq. (|58| ) into Eq. ( p9[) then gives 



Eii(q.w) = J—, rEn(q, 



bw) (60) 



and, accordingly, 



s ctt (q^) 



'2„ 



z D (q, w) 



(61) 



frc- 



where Z n (q, w) = l/g n (q, w) is momentum and 
qucncy dependent sheet impedance. 

The form of Eq. ( |6l"| ) suggests that the response of a 
current in a two-dimensional system to external electric 
field is the same as if the proper sheet impedance of a 
system Z a (q, <jj) has been connected in series with the 
other contribution, which can be considered as the ef- 
fective impedance of the empty space surrounding this 
system. This additional contribution is purely inductive 
and corresponds to 



L 8 (q) 



Mo 
2q 



(62) 



In the case of a superconducting system in the low 
frequency limit Z n (q, w) w —iu>L n , where L a is the ef- 
fective sheet inductance, substitution of which into Eq. 
( |60| ) allows to rewrite it as 



E ||(q,w) 



Aq 



1 + Aq 



En(q,w) 



where 



A 



2Lg 

Mo 



(63) 



(64) 



is the (two-dimensional) magnetic field penetration 
length |2|] already discussed in Sec. 3. Since in the con- 
sidered system electric field appears only due to presence 
of ac magnetic field, the same length describes as well the 
screening of the electric field. 

The form of the fluctuation-dissipation theorem ob- 
tained after substitution of Eq. ( |6l| ) into Eq. ( p^ ) 



i'i'Jqu = 2TRe- 



1 



-iuL s (q) + Z a (q,uj) 



(65) 



being completely independent of the details of the struc- 
ture of the particular system used for its derivation, one 
can expect it also to be valid for other two-dimensional 
superconducting (or simply conducting) systems, in par- 
ticular superconducting films. 



VI. RESULTS AND DISCUSSION 

Substitution of Eq. ( |65| ) into Eq. ( p^ ) gives the expres- 
sion for the flux noise power spectrum 



HIT f d 2 q 



(2tt) 



2 F 2 (q)Re 



1 



-iu>L s {q) + Z n (q,u>) 



(66) 



which is the central result of this work. It allows instead 
of constructing special theories explaining frequency de- 
pendence of S(u>) in different regimes to relate it with 
already known properties of Z u (q,u>). 

It is of interest to compare Eq. (^) with the expres- 
sion for the quantity [the correction to frequency depen- 
dent mutual impedance Z m (u)], which is measured in the 
framework of the two coil method and defined as the ra- 
tio of the measured voltage to the driving current. In our 
notation this expression JlOj] can be rewritten as 



SZ m {ui) = -u 2 4 



F{q)F{q) 



(2n) 2 -iu;L s (q) + Z n (q,co) 



(67) 



where F(q) is the geometrical factor analogous to F(q), 
but describing the driving coil. Note that applicability 
of Eq. ( |66[ ) and Eq. ( |67j ) for description of a Josephson 
junction array requires the lattice constant of the array 
to be much smaller then the geometrical parameters of 
the coil(s) (r and h) entering the expressions for F(q) 
[and %)]. 

Comparison of Eq. ( |67| ) with Eq. ( |66| ) shows that the 
real part of 8Z m (io) is determined by the same (real) com- 
ponent of g cS (q,uj) as the noise spectrum, so in absence 
of a difference between F(q) and F(q) the noise spectrum 
S(u>) would be simply proportional to Re[6Z m (tjj)]: 



2T 

S(cj) j Re[SZ m ( 



(68) 



In the case of the simplest coil (a single circular loop 
of radius r) Eq. (pq) is reduced to 



J 2 (rq) 



S(u>) = TT^Tr I dq 

Jo 11 



x Re 



exp(-2hq) (69) 
1 



-iu>L s (q) + Z a (q,u>) 



The analogous equation derived by Kim and Minnhagen 
P7[ differs from Eq. ([39|) basically (i) by the absence of 
the term L s (q) and (ii) by the presence in the integrand 
of the additional factor q 4 . The former of the two discrep- 
ancies is rather natural, since in Ref. jlTj the screening 
effects have not been taken into account, whereas the 
latter we believe to be the consequence of the incorrect 
calculation of the magnetic field produced by currents in 
the array. 



6 



In Ref. |Tij this magnetic field is calculated as the sum 
of magnetic fields produced by current loops associated 
with lattice plaquettes, whereas the magnitude of a cur- 
rent in each loop is assumed to be given the directed sum 
of the currents in the junctions 7 na along the perimeter 
of a plaquette (in the present work this sum is denoted 
V x 7). In such procedure the current of each junction 
is counted twice (as giving contributions to the loop cur- 
rents associated with the two neighboring lattice plaque- 
ttes) but with opposite signs, so these two contributions 
almost cancel each other, which leads to appearance of 
the additional (/-dependent factor in comparison with our 
Eq. (69). In a more consistent implementation of this ap- 
proach the values of loop currents associated with arrays 
plaquettes should be chosen in such a way that the value 
of the current on each junction is given by the difference 
of the loop currents associated with the two neighboring 
plaquettes. These so called mesh currents |2l]]24j 7™ are 
related with the currents on the junctions I na as 



7- = -A£ X (V x 7) r 



(70) 



which explains the appearance of the additional factor of 
q 4 in calculation which uses (V x 7) n instead of 7™. 

The main contribution to the integral in Eq. (B9) is 
coming from the region 



< q < 1/h, 



(71) 



so in the situation when h exceeds all the microscopic 
scales responsible for the q dependence of Z a (q,u>), one 
can replace in Eq. ( p9|) Zn(q,u>) by 

Z D (w) = —iuiL a (lo) + R n (w) = lim Z a (q, w) (72) 

and use S(u>) to extract information about L a (to) and 
Rn(ui). 

Let us first discuss the limit when the effects of screen- 
ing can be neglected. This is possible if in the essential 
part of the interval ( |7l| ) one can neglect L s (q) in com- 
parison with Lo (w) and therefore requires h <C A (and 
not ft> Aas has been claimed in Ref. 0] and Ref. [jL7| ) . 
In that case Eq. (p9) is reduced to 



S(w)=iri4Tr 2 Y(h/r) 



where 



Y(h/r) 



4irq) 



dq 

q 

l 

2 



exp (— 2hq) 



(73) 



(74) 



for /i/r<l 
^(£) 2 fartyr>l ■ 

This means that with increase of h the amplitude of the 
noise has first to change very slowly and then to decay 



proportionally to 1/h 2 . The experimental results of Fes- 
tin et al. obtained on thin YBCO film at h/r ~ 1, are 
compatible with 1/h 2 dependence even better than with 
1/h 3 dependence predicted by Kim and Minnhagen [fj"7| . 

In the opposite limit of strong screening (A <C h) 
Lo(u>) is negligibly small in comparison with L s (lo) and 
S(uj) depends only on R n (co), but not on Ln(u>). In par- 
ticular, for h <^ir integration in Eq. (|69|) gives 



S(w) 



It Bo (oQT 
h ZP 



for lo <C 
for lo 3> 



Hoh 
Mop 
Hoh 



(75) 



Note that for A <c h <C r and small lo one obtains the 
dependence S(lo) oc 1/lo, which is observed in noise mea- 
surements in Josephson junction arrays [[^] over four 
decades in frequency. 

The idea that this particular dependence appears be- 
cause SQIUD integrates contributions from wide interval 
of wavevectors has been put forward by Wagenblast and 
Fazio HU . Our analysis has shown that in order to obtain 
the 1/lo dependence over four decades in frequency one 
should have r/h > 10 4 , which in the experiments of Refs. 
j?],H definitely has not been fulfilled. Thus the origin of 
the 1 /lo dependence of the flux noise power spectrum ob- 
served experimentally in Josephson junction arrays (in 
particular, at tempratures lower then the expected tem- 
perature of the Berezinskii-Kosterlitz-Thouless (BKT) 
phase transition p5| , p6[ |) still remains to be elucidated. 

Substitution of the results of Ambegaokar et al. Gffl 
for vortex pair contribution to impedance into Eq. (|66| ) 
or its simplified analogs allows to reproduce the depen- 
dence S(lo) oc 1/lo only exactly at the temperature of 
the BKT transition. On the other hand, the analogous 
frequency dependence of the noise spectrum observed in 
thin YBCO films || is ascribed to vortex hopping be- 
tween neighboring pinning centers, a finite vortex con- 
centration being associated with presence of a residual 
magnetic field. The same mechanism may be responsible 
for the results of experiments on Josephson junction 
arrays. 
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